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Thursday,  December  12,  1968 

4:00-  8:00  P.M.  -  Registration 

Lobby,  Center  for  Continuing  Fducation 


Friday,  December  13,  1968 
8:90-  9:00  A.M.  -  Registration 

Lobby,  Center  for  Continuing  Education 

9:00-10:10  ?,*»,  -  General  Session  1 

Auditorium,  Center  for  Continuing 
Education 

Chairman:  Dr.  Charles  E.  Carver,  Jr. 
Co-Manager,  Project  THESIS, 
University  of  Massachusetts,  Amherst 


9:00-  9:15  A.M. 


9:15  A.m. 


10:10-10:30  A.**. 
10:30-12:00  P.m. 


10:30-11:00  A.M, 


11:00-11:30  A.M. 


-  Welcome:  Dr.  Frederick  D.  Rossini, 

Vice-President  for  Pesearch  and. 

Special  Programs, 

University  of  Notre  Dame 

-  Lecturer  ,*DEEP  SUBMERGENCE  VEHICLE 

DEVELOPMENT*'"/  ' 

*’r.  Roger  D.  Fuller,  Program  Manager, 
Deep  Submergence  Search  Vehicle, 
Lockheed  ’’issiles  and  Space  Company, 
Sunnyvale,  California  / 

/ 

-  Coffee  Break  -  Patio 

-  Technical  Session  fA  &  IB  (Concurrent) 

Technical  Session  IA 
Room  112,  Center  for  Continuing 
Education 

Chairman:  Dr.  S.  M.  Bember., 

Department  of  Civil  Engineering 
University  jf  Massachusetts 

-^STATISTICAL  ERROR  ANALYSES  NAVIGATION 
SYSTEMS 

C.  E.  Hutchinson, 

_ University  of  .Massachusetts 

-^^Sensitivity  of  stability  of  a  class  of 

NON-LINEAR  TI.MEVARYING  NETWORKS 
R.  W.  Liu,  University  of  Notre  Dame  "3 
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Friday ,December  13,  1963  (cont'd.) 
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11:30-12:00  P.*<.  -  /"SUBSYSTEM  VALUE  AND  COST-DEEP  OCEAN 

SYSTEMS*'  / , 

W.  E.  Heronelnus, 

University  of  Massachusetts 

Technical  Session  IB 

Room  210,  Center  for  Continuing  Education 
Chairman:  Dr.  K.  T.  Yang, 

Departr.ent  of  Mechanical  Engineering 
University  of  Notre  Dame 

10:30-11:00  A.M.  -  H’URBULEMT  FREE  CONVECTION*^ 

A.  A.  Szewczyk,  v 

University  of  Notre  Dame  "'" 

11:00-11:30  A .  M  .“^SURFACE  WAVE  ATTENUATION  FOR  OCEAN 

ENGINEERING  PURPOSES/'"  J 
J..JI.  Cornell,  University  of  Massachusetts 

11:30-12:00  P.M.  ^  /DRAG  EFFECTS  OF  POLYMER  ADDITIVES  IN 

TURBULENT  NON-NEWTONIAN  FLAT  PLATE 
BOUNDARY  LAYER  FLOW**''"  /  \ 

R.  L.  Sampson,  University  of  Notre  Dame 

12:15-  2:00  P.M.  -  Luncheon,  Center  for  Continuing  Education 

Presiding:  William  B.  Walker,’  Commander , 

USN,^" 

Director  Undersea  Programs, 

Office  of  Naval  Pesearch,  Washington,  D.C. 


Address:  >THE  AIMS  AND  OBJECTIVES  OF 
PPOJFCT  THEMIS^  -  - 


PPOJFCT  THEMISr  -  ( 

Dr.  Arwin  A.  Dougal,  Assistant 
Director  (Research)  and  Chairman, 
Project  THEMIS  Task  Group,  Office  of 
the  Director  of  Defense  Research  and 
Engineering,  Washington,  L>.  C. 

2:00-  4:00  P.M.  -  Technical  Session  IIA  &  IIB  (Concurrent) 

Technical  Session  IIA 

Poom  112,  Center  for  Continuing  Education 
Chairman:  Dr.  C.  E.  Hutchinson, 

Department  of  Electrical  Engineering 
University  of  Massachusetts 
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2:00- 

2:30- 

3:00- 

3 : 30- 

4:00- 

2:00- 

2:30- 

3:00- 

3:30- 


2 • 30  P.M.  -  <*THE  VERTICAL  HOLDING  CAPACITY  OF  MARINE 
ANCHORS  IN  SANDS/y 
S.  PI.  Bemben  * 

_  University  of  Massachusetts 

3-00  P.M.  -  /FATIGUE  STRENGTH  OF  SODA-LIP 'E-SILICA 
GLASS^V) 

J.  E.  Ritter,  Jr. 

University  of  Massachusetts 

3:30  P.M. /THE  EFFECT  OF  STRESS  ON  MAGNETOSTRICTION  ^ 
C.  W.  Allen 

University  of  Notre  Dame 

4:00  P.^.  ^/DEVELOPMENT  OF  HIGH  MAGNETOSTRICTION 

TRANSDUCER  ALLOYS^-  • _ 

A.  E.  Miller  '  \ 

University  of  Notre  Dame 

4:15  P.M.  -  Coffee  Break  -  Lobby 
Technical  Session  IIB 

Room  210,  Center  for  Continuing  Education 
Chairman;  Dr.  A.  A.  Szewczyk 
Department  of  Mechanical  Engineering 
University  of  Notre  Dame 


2*30  P.M.  -  /6n  THE  FLOW  FIELD  APOUND  A  MARINE 
DUCTED  PROPELLOW^y 
D.  E.  Cromack, 

University  of  Massachusetts 

3:00  P.M.  ^  ^PHYSICAL  AND  NUMERICAL  EXPERIMENTS  IN 
HYDRODYNAMIC  WAKE  FLOWS /•"/ 

T.  J.  Mueller  J  \ 

University  of  Notre  Dame 

3-.30  p.M.  -  /Transmission  of  sound  through  wipe 
SCREENS  IN  A  STEADY  FLOW®-''/ 

G.  A.  Russell  y 

University  of  Plassachusetts'''^ 

A  ;00  P.M.  -  /tAMINAR  FREE  CONVECTION  FROM  A  NON- 
ISOTHERMAL  VERTICAL  PLATE  IN  A 
TEP'PERATURE-STPATIFIED  ENVIRONMENT*/ 
K.  T.  Yana 

University  of  Notre  Dame 


4:00-  4:15  P.M.  -  Coffee  Break  -  Lobby 


4:15-  5*10  P.M. 


5*45-  6:^5  P.M. 


G  :  45-  C : 30  P.M. 


P : 30  P.M. 


General  Session  II 

Chairman:  Professor  W.  E.  Heronenus 

Department  of  Civil  Enqineerinq 
University  of  Massachusetts 
Lecture:  ^RESEARCH  NEEDS  FOR  THE  DEEP 
OCEANS'''/  v 
Dr.  Allyn  C.  Vine, 

Senior  Oceanographer, 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts 

Coe}* tail  Reception, 

Mahogany  Room  -  Morris  Inn 

Symposium  Banquet/  "’orris  Inn 
Presiding:  Dr.  Joseph  C.  Hogan 
Dean,  College  of  Engineering 
University  of  Notre  Dame 

General  Session  III 

Auditorium,  Center  for  Continuing  Education 
Chairman:  Dr.  Harry  C.  Saxe 

Program  Manager,  Project  THEMIS 
University  of  Notre  Dame 
Lecture  f*  /"CONSTRUCTION  AND  PROPOSED 
USFS  OF  THE  SUBMERSIBLE- 
BEN  FRANKLIN?"'  /  , 

Walter  R.  Scott,  Jr. 

Director,  Ocean  Systems  Department 
Grumman  Aircraft  Engineering  Corporation 
Bethpage,  New  York. 


aturdey,  December  14,  1968 


9:09-  9:55  A.D.  -  General  Session  IV 

Auditorium,  Center  for  Continuing  Education 
Chairman:  Dr.  L.  H.  N.  Lee 

Department  of  Engineering  Science 
University  of  Notre  Dame 
Lecture:*  *t>EFPSTAR  4000  -  AN  OCEANO¬ 
GRAPHIC  RESEARCH  TOOL*"/  ^ 

*’r.  Richard  Usry  „*■ 

Deeps tar  Marketing  Representative 
Westinnhouse  Underseas  Division 
Ocean  Research  and  Engineering  Center 
Ealtimcre,  Maryland 


9:55-10:15  A.M.  -  Coffee  Break  -  Patio 


Saturday,  December  14,  1968  (cont'd.) 

10-15-11:45  A.M.  -  Technical  Session  III?  &  IIIB  (Concurrent) 

Technical  Session  IIIA 

Pcom  112,  Center  for  Continuing  Education 
Chairman:  Dr.  B.  D.  Cullity 

Department,  of  ^etalluraical  F.ngi  nociiny 
and  "aterials  Science 
University  of  Nctre  Dame 

• 

10:3  5-10:^5  A  .M.  -  “THE  APPLICATION  OF  DYNAMIC  RELAXATION 

TO  THE  FINITE  ELEMENT  METHOD  OF 
STRUCTURAL  ANALYSIS'*  • 

S.  Kelsey 

University  of  Notre  Dame 

10-45-11-15  a.m.  -^ton-linear  analysis  of  platfs  and  shflls" 

M.  M.  Miller 

University  of  Massachusetts 

^  a 

11:15-11:45  A.M.  -  -•‘INELASTIC  AXISY*'MFTRIC  BUCKLING  OF  RING 

STIFFENED  CYLINDRICAL  SHELLS'* 

L.  H.  N.  Lee 

University  of  Notre  Dame 
Technical  Session  IIIB 

Room  210,  Center  for  Continuing  Education 
Chairman:  Dr.  J.  F.  Ditter,  Jr. 

Department  of  Mechanical  and  Aero- 
Space  Engineerinc 
University  of  Massachusetts 

■V 

10:15-10:45  A.M.  -  “LICET  SCATTERING  AS  METHOD  OF 

DETERMINING  THE  ACOUSTICAL  AND  TRANSPORT 

PROPERTIES  OF  FLUIDS* 

•  ^  1 
A.  A.  Morhewicz 

University  of  Notre  Dame 
\  , 

10  4  5-11:15  A.f*.  ‘"'DYNAMIC  P^OGFAThllNO  AND  THE  DISCPETE 

MAXIMUM  PRINCIPLE*  J 
P.  J.  LeaJ-e 

University  of  Notre  Dame 

11:15-11:45  A.M.  ^>NlCPOSTRUCTUnE  EFFECTS  IN  THE  ANALYSIS 

OF  COMPOSITES®"",’  ^ 

T.  Ariman 

University  of  Notre  Dame  " 


12:00-  1.45  P.N. 


1:45- 


-  Luncheon 

Center  for  Continuing  Education 
Presidinq:  Dr.  C.  E.  Carver,  Jr. 

Co-Manager,  Project  THEMIS 
University  of  Massachusetts 
Lecture  v’  "SEALAE  -;HErf 
John  P.  Lindsay,  Caotain  USN 
Deputy  Project  .Manager 
Deep  Submergence  Systems  riojeot- 
Chevy  Chase,  Maryland 

3:30  p.M.  -  Closing  Session 

Auditorium,  for  Continuing --Education 
Panel  Discussion  -  Tonic:  "WHAT'S  AHEAD 
IM  OCEAN  ENGINEERING"  ,  . 

Moderator:  Dr.  Kenneth  G.  Picha 

Program  Manager,  Project  THEMIS 
and  Dean,  School  of  Engineering 
University  of  Massachusetts 

Panel  Consisting  of: 

Dr.  Victor  C.  Anderson 
Acting  Director,  Marine  Physical 
laboratory 

ScriDps  Institution  of  Oceanography 
La  Jolla,  California 

Dr.  Frank  A.  Andrews 

Director,  Institute  of  Ocean  Engineering 
The  Catholic  University  of  America 
Washington,  D.C. 

Captain  Henry  a.  Arnold 
Marine  Sciences  Council 
Executive  Office  of  the  President 
Washington,  D.C. 

Commander  William  B.  Walker 
Director,  Undersea  Proarams 
Office  of  Naval  Research 
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STATISTICAL  ERROR  ANALYSES  OF 
NAVIGATION  SYSTEMS 

bv 

Charles  F.  Hutchinson 
Electrical  Encrineering  Department 
University  of  Massachusetts 
Amherst,  Massachusetts 

The  problem  of  accurate  navigation  is  critical  to  any 
marine  operation.  Specifically  one  would  like  to  have  precise 
continuous  position  information  in  all  tvoes  of  weather  at  all 
operating  locations.  This  of  course,  is  not  possible,  but  it 
is  a  commendable  goal  for  any  navigation  system.  Due  to  the 
importance  of  navigation  and  the  expense  of  obtaining  precise 
position  information  considerable  effort  is  expanded  trying  to 
analyze  the  errors  inherent  to  any  svstem.  Fistoricallv ,  how¬ 
ever,  it  is  only  recently  that  the  problem  of  error  analyses 
of  navigation  systems  has  been  approached  from  a  statistical 
viewpoint. 

This  paper  presents  the  concent  of  error  analysis,  both 
deterministic  and  stochastic.  The  procedures  are  outlined 
using  the  inertial  system  as  an  example,  however,  the  tech- 
nigues  have  application  to  anv  navigation  system. 

The  task  of  analyzing  errors  in  navigation  systems  may  be 
illustrated  by  the  following  diagram. 


That  is,  three  separate  sub-problems  exist.  Starting  from  the 
right  in  the  diagram  above,  the  first  task  is  to  determine 
what  svstem  outputs  are  critical  to  the  system  mission,  For 
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most  applications  this  would  simply  be  latitude  and  longitude. 
Although  there  are  cases  where  other  outputs  become  important 
as  well.  For  the  sake  of  example  we  will  consider  latitude 
error  in  an  inertial  system,  as  the  specific  output  of  concern. 

The  next  phase  of  the  analysis  is  to  determine  the  sources 
of  error  and  then  the  dynamics  which  describes  how  these 
sources  effect  the  outputs.  For  any  system  of  practical  value 
there  will  be  many  such  error  sources.  ^or  the  inertail  sys¬ 
tem  the  list  might  include  such  items  as: 

gyro  drift  rates 
accelerometer  errors 
velocity  reference  errors 
platform  misalignments 
initialization  errors 
etc . 

Each  of  these  error  sources  might  in  turn  be  further  divided. 

Such  as,  under  gyro  errors  we  might  have: 

mass  unbalance 
temperature  effect 
g  sensitive  drift 
etc. 

The  dynamics  involved  are  normally  determined  by  noting 
the  physical  properties  of  the  system.  The  various  components 
must  be  analyzed  to  determine  their  individual  dynamic  behavior. 
This  portion  of  the  problem  is  probably  the  most  readily  avail¬ 
able  to  the  analyst,  '"'uch  work  has  been  done  in  this  area. 


'T~«r  - 
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SENSITIVITY  OF  STABILITY  OF  A  CLASS  OF 
NON-LINEAR  TIMEVARYING  NETWORKS 

by 

R.  W.  Liu 

Electrical  Engineering  Department 
University  of  Notre  Dame 
Notre  Dame,  Indiana  46556 

Sensitivity  problems  of  stabilities  of  network  models  are 
considered.  There  are  two  reasons  to  choose  a  nonlinear  net¬ 
work  model,  instead  of  a  nonlinear  systems  model  (state  equa¬ 
tions)  .  Firstly,  since  any  nonlinear  system  equation  can  be 
synthesized  by  a  nonlinear  network,  no  generality  is  lost  by 
studying  nonlinear  networks.  Secondly,  unlike  system  models, 
a  network  has  a  structure:  interconnections  of  components. 
Therefore,  it  is  more  convenient  to  consider  sensitivity  pro¬ 
blems  in  a  network  model. 

The  problem  is  to  study  the  effect  on  the  stability  behav¬ 
ior  of  a  network  with  a  variation  of  a  (nonlinear)  component. 
Two  kinds  of  stabilities  are  studied:  the  global  asymptotic 
stability  with  zero-input  and  the  convergence  of  solutions 
with  an  input. 

Two  classes  of  networks  are  considered:  passive  net¬ 
works  and  two-element-kind  networks.  Necessary  and  sufficient 
conditions  for  each  case  are  obtained.  The  conditions  are  in 
terms  of  network  topology  and  characteristic  equations  of 
network  elements.  Therefore,  it  is  very  easy  to  test. 
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SUBSYSTEM  VALUE  AND  COST-DEEP  OCEAN 
SYSTEMS 

by 

William  E.  Heronemus 
Civil  Engineering  Department 
University  of  Massachusetts 
Amherst,  Massachusetts 

This  paper  is  intended  to  describe  the  total  scope  of  a 
University  of  Massachusetts  THEMIS  project  in  which  four  princi¬ 
pal  investigators  and  five  graduate  assistants  are  conducting 
research  as  individuals,  but  all  correlating  with  one  central 
project  theme.  The  project  title  in  itself  is  deliberately 
broad:  "Deep  Submergence  Vehicle  System  Design."  We  are  in¬ 

terested  in  the  DSV  System  in  its  broadest  and  narrowest  con¬ 
notations.  We  are  interested  in  the  designs  of  others,  and  we 
fully  intend  to  create  many  designs  of  our  own.  We  are  con¬ 
ducting  research  in  the  highly  specialized  field  of  subsystems 
and  components  for  navigation,  guidance  and  control,  but  we 
want  to  evaluate  them  and  the  value  of  our  research  only  in 
terms  of  what  they  mean  to  the  cost  and  value  of  a  Deep  Ocean 
System. 

This  is  the  first  attempt  at  a  summarizing  report  on  the 
entire  research  project.  It  is  a  very  early  report:  the  pro¬ 
ject  has  scarcely  been  started  from  the  broad  point  of  view 
albeit  it  is  a  year-and-a-half  alono  in  some  of  its  specialized 
efforts  to  date  are  reported  h**  Hutchinson  and  Monopoli  et  al. 
in  several  other  THEMIS  repo,.  one  of  which  will  be  read 
here  at  this  meeting. 
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The  Deep  Ocean  System  is  thus  a  working  system,  a  mission 
and  task  accomplishing  system,  not  an  inanimate  object.  Seven 
examples  of  major  working  devices  which  could  be  found  in  the 
02  tier  for  a  Deep  Ocean  System  are: 

(a)  A  Deep  Submergence  Vehicle — any  one  of  many  different 
varieties,  depending  upon  the  Mission  and  Tasks  of 
the  Deep  Ocean  System 

(b)  An  Off-shore  Oil  Drilling  Rig 

(c)  An  Off-shore  Oil  Production,  Storage  and  Distribution 
Plant 

(d)  A  Fish  Food  Processing  Ship 

(e)  A  Sea  Bed  Mining  Dredge 

(f)  An  Observational  Laboratory  for  Marine  Sciences 
Education 

(g)  A  "Sea  Lab" 

The  current  University  of  Massachusetts  THEMIS  Project  limits 
the  02  tier  to  a  Deep  Sea  Vehicle. 

Results  to  Date 

This  project  in  its  broad  sense  has  just  been  started. 

•  Results  to  date  are: 

(a)  Verification  that  the  one,  t**o  and  three  sphere 
"hulls"  are  best  suited  for  the  initial  scope  of  the 
investigation. 

(b)  Determination  of  a  set  of  sphere  diameters  appropriate 
to  one  through  five  man  occupancy. 

(c)  Calculations  that  demonstrate  that  the  three  selected 
energy  storage  subsystems  are  reasonably  competitive. 


IA  (6) 


(d)  Verification  that  the  project  should  neither  ignore 
the  lower  depth  bands  nor  the  lower  cost  less  so¬ 
phisticated  pressure  hulls. 

(e)  Verification  that  the  many  published  computer  pro¬ 
grams  for  naval  architectural  computations  are  not 
directly  applicable  to  computation  of  the  DSV  hull 
which  is  a  chain  of  displacement  shaoes  enveloped  by 
a  free-flooding  outer  hull. 

(f)  Collection  and  partial  analysis  of  considerable  data 
describing  existing  DSV's,  their  subsystems  and  their 
components. 

(g)  Initial  computations  with  a  DSV  Simulation  Computer 
Program. 

(h)  Formulation  of  policy  concepts  for  the  kinds  of  navi¬ 
gation,  command  and  control  subsystems  to  which 
emphasis  will  be  applied. 

(i)  Initial  description  of  kinds  of  variable  ballast  and 
trim  subsystems  to  which  emphasis  will  be  applied. 

Plans  for  Future  Work 

During  the  spring  and  summer  considerable  effort  will  be 
expended  on  the  mathematical  and  computer  aids  required  on  the 
project.  The  collection  and  analysis  of  data  defining  existing 
DSV's  will  continue.  The  collection  and  analysis  of  subsystem 
and  component  data  and  synthesis  of  subsystems  for  the  initial 
families  of  DSV's  will  continue.  Synthesis  and  graphic  delin¬ 
eation  of  the  initial  families  of  DSV's  will  be  started.  De- 
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finition  of  subsystem  value  and  cost  parameters  will  continue. 

By  iteration  those  subsystems  found  to  be  of  most  value  will 
become  ever  more  apparent.  The  possibilities  for  more  versa¬ 
tile  Deep  Sea  Vehicles  of  lower  initial  and  ongoing  cost  will 
follow  the  identification  of  value  and  cost  of  their  subsystems. 

The  long-range  goal  is  that  any  new  subsystem,  component, 
material,  fabrication  process,  any  research  product,  aimed 
toward  improved  Deep  Submergence  Vehicles,  can  be  evaluated 
for  value  and  cost  in  any  of  a  very  large  number  of  DSV  fami¬ 
lies.  The  long-range  goal  is  a  continuing  trade-off  study  in 
the  entire  field  of  DSV  technology. 
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TURBULENT  FREE  CONWCTION 
by 

A.  A.  Szewczyk 

Denartment  of  Mechanical  Engineering 
University  of  Notre  Dame 
L  Notre  Dame,  Indiana 

The  engineering  relevance  of  turbulent  free  convection 
on  heated  surfaces  in  connection  with  under-sea  heat  rejection 
systems  is  first  emphasized.  The  talk  describes  a  different* a i 
phenomenological  theory  for  turbulent  free  convection  along 
a  heated  vertical  plate.  .  The  total  viscosity 

is  taken  to  be  governed  by  a  rate  equation  which  takes  into 
account  turbulent  convection,  diffusion;  generation  and  decay, 
while  the  turbulent  thermal  diffusivity  is  based  on  a  constant 
turbulent  Prandtl  number.  Calculations  based  on  experimental 
data  for  air  show  that  the  turbulent  boundary  layer  consists 
of  a  viscous  sublayer  with  nonlinear  velocity  distribution, 
an  inner  layer  dominated  by  turbulent  generation  and  decay 
and  sensitive  to  the  Grashof  number,  and.  an  outer  layer 
governed  by  turbulent  convection  and  diffusion,  and  that 
Corrsin's  viscous  super  layer  also  exists  in  the  nresent  case. 

Also  described  is  the  progress  on  the  accompanying 
experimental  study  in  which  temporal  and  spectral  details  of 
the  turbulent  boundary  layer  of  air  on  a  cooled  vertical 
plate  23  feet  and  6  inches  in  height,  are  to  be  measured. 

The  primary  purpose  here  is  to  provide  much  needed  data  to 
further  substantiate  the  theoretical  findings . 


f 
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SUPPACE  WAVE  ATTENUATION  FOR 
OCEAN  ENGINEERING  PURPOSES 

by 

Dr.  C.  E.  Carver,  Jr. 

Dr.  J.  M.  Colonell 
Department  of  Civil  Engineering 
University  of  Massachusetts 
?mherst,  "assachusetts 

The  lauching  and  recovery  phases  of  submersible  vehicle 
operation  can  be  extremely  hazardous  to  personnel  and  the 
vehicle  itself  when  surface  vrave  activity  exceeds  specified 
operating  limits.  To  minimize  dangers  produced  by  high  seas 
and  the  resulting  excessive  vehicle  motions,  it  is  suggested 
that  mobile  breakwater  devices  be  deployed  by  the  mother 
ship.  The  primary  objective  of  this  investigation  is  to 
evaluate  the  feasibility  of  using  either  hydraulic  or  pneumatic 
breakwater  devices  to  create  an  acceptable  environment  for 
the  lauching  and  recovery  of  a  submersible  vehicle. 

Previous  studies  of  these  breakwater  devices  are 
summarized  and  a  brief  of  the  associated  hydrodynamical  theory 
is  presented.  The  experimental  orocrram  for  this  investigation 
will  be  conducted  in  a  wind  and  water  wave  research  facility 
now  under  construction.  This  facility,  which  has  a  cross- 
section  that  is  2  feet  wide  by  4  feet  high  and  has  a  length 
of  48  feet,  will  enable  the  study  of  waves  generated  by  winds 
with  speeds  up  to  60  fps.  The  experimental  schedule  for  the 
investigation  is  discussed  with  particular  attention  to  the 
proaram  objectives. 
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DRAG  EFFECTS  OF  POLYMER  ADDITIVES 
IN  TURBULENT  MON-NEWTONIAN 
FLAT  PLATE  BOUNDARY  LAYER  FLOW 

by 

R.  L.  Sampson 

Department  of  Chemical  Engineering 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

The  current  status  on  drag  reduction  with  polymer  addi¬ 
tives  in  pine  flows  as  well  as  in  boundary  layer  flows  is 
briefly  reviewed.  This  THEMIS  project  deals  essentially  with 
the  experimental  study  of  drag  reduction  effects  of  polymer 
additives  in  the  turbulent  boundary  layer  on  a  flat  plate. 

A  fully  instrumented  water-tunnel  facility,  specially  designed 
for  this  experimental  study  with  suitable  control  in  flow 
and  temperature,  has  been  constructed  and  will  become 
operational  shortly.  A  laser  velocimeter  system  is  developed 
and  its  measurements  have  been  compared  with  that  of  a  pitot 
tube. 

A  predictive  analysis  of  drag  coefficients  for  turbulent 
non-Newtonian  flat  plate  boundary  layer  flow  has  also  been 
developed.  With  computer  assistance,  flat  plate  drag  co¬ 
efficients  are  calculated  as  functions  of  a  generalized 
Reynolds  number  and  fluid  rheological  properties.  In  utilizing 
both  the  laser  velocimeter  and  pitot  tubes  for  velocity 
measurements  and  the  capillary  tube  viscometer  for  deter¬ 
mination  of  rheological  properties,  data  for  turbulent 
boundary  layer  flow  will  be  obtained  for  water  and  aqueous 
carbopol  solutions.  Dilute  carbopol  solutions  are  char¬ 
acterized  by  non-Newtonian  power  low  behavior  with  negligible 
elasticity.  These  experimental  data  will  be  compared  with 
the  theoretical  predictions. 
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THE  VERTICAL  HOLDING  CAPACITY 
OF  MARINE  ANCHORS  IN  SAND 

by 

S.  M.  Bemben 

Department  of  Civil  Engineering 
University  of  Massachusetts 
Amherst/  Massachusetts 

A  large-scale  outdoor  test  bin  facility  has  been 
constructed.  The  appurtenances  include  a  well  point  system 
for  drawdown  and  recharge  purposes  and  a  movable  reaction 
frame.  Pulling  loads  are  applied  by  a  hydraulic  jack  and 
chain  assembly  and  are  measured  by  strain  gage  load  cells. 

To  date,  model  anchors  having  projected  horizontal  areas 
of  about  1.5,  15,  30,  60  and.  120  square  inches  have  been 
inserted  to  various  depths  ranging  between  2  feet  and  6  feet 
in  a  flooded  sand  soil  and  load  tested  for  vertical  holding 
capacities.  For  insertion,  the  anchors  are  attached  to  the 
lower  end  of  a  rigid  pipe  and  vibrated,  into  place  by  means 
of  a  vibrator  unit  attached  to  the  upper  end.  A  horizontally 
cancelling  double  eccentric  vibrator  and  a  single  eccentric 
vibrator  have  been  employed;  both  units  are  driven  by  a 
controlled  variable  speed  motor. 

Preliminary  test  data  indicate  that,  for  the  anchors  test¬ 
ed  in  the  one  sand  used,  an  optimum  resisting  load  develops 
when  the  ratio  of  the  depth  of  burial  to  the  square  root  of 
the  horizontal  projected  area  of  the  anchor  is  about  5  to  1. 

The  ratio  of  the  peak  resisting  loads  to  horizontal  projected 
area  appears  to  increase  slightly  with  increasing  anchor  size. 
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FATIGUE  STRENGTH  OF  SODA-LIME-SILICA  GT.ASS 

by 

J.  E.  Ritter,  Jr. 

Aerospace  &  Mechanical 
Engineering  Department 
University  of  Massachusetts 
Amherst,  Massachusetts 

The  effect  of  strain  rate  on  the  room-temperature  strength 
of  soda-lime-silica  glass  was  determined.  The  strength  of 
acid-etched  glass  rods  increased  continuously  from  190  200 
to  284  800  psi  on  increasing  the  strain  rate  from  9X10  4  to 
15X10  2in./in./min.  No  difference  in  strengths  was  observed 
on  testing  in  laboratory  atmosphere  or  with  the  specimens 
wetted  with  distilled  water.  For  abraded  glass  rods  the 
strength  in  laboratory  atmosphere  varied  from  12  555  to  16  540 
psi  on  varying  the  strain  rate  from  3X10_3  to  9X10  2 in. /in./ 
min.  Tests  conducted  with  the  abraded  samples  wetted  with 
distilled  water  caused  approximately  a  10%  reduction  in 
strength,  but  no  appreciable  change  in  the  strain-rate 
sensitivity.  The  stress  corrosion  model  proposed  by  R.  J. 
Charles  successfully  predicts  the  room  temperature  strain- 
rate  sensitivity  of  the  failure  process  for  both  acid-etched 
and  abraded  glass  samples.  Using  relationships  derived  from 
this  model,  it  was  further  shown  that  the  stress  concen¬ 
tration  relationship  proposed  by  Inglis  is  valid  for  acid-etched 
glass . 


. 
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THE  EFFECT  OF  STRESS  ON  MAGNETOSTRICTION 

by 

C.  W.  Allen 

Metallurgical  and  Materials  Science  Engineering 
University  of  Notre  Dame 
Notre  Dame,  Indiana  46556 

An  elementary  introduction  to  a  number  of  facts  regarding 
magnetostriction,  including  its  anisotropy,  is  presented.  The 
effect  of  stress  on  magnetostriction  enters  in  terms  of  the 
effect  of  stress  on  the  distribution  of  magnetization,  a  change 
of  which  in  turn  modifies  the  magnetostriction.  This  is  illus¬ 
trated  by  experimental  curves  showing  the  effect  of  apolied 
tensile  and  compressive  loads  on  the  magnetization  curves  for 
iron  and  nickel  of  commercial  purity  and  on  the  technical 
magnetostriction  curves  for  these  situations.  Preliminary  data 
are  presented  on  magnetostriction  of  iron  and  nickel  with  very 
low  frequency  magnetic  field  cycling  at  several  levels  of 
applied  stress  and  average  magnetization. 
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DEVELOPMENT  OF  HIGH  MAGNETOSTRICTION 
TRANSDUCER  ALLOYS 

by 

A.  E.  Miller 

Metallurgical  &  Materials  Science 
Engineering  Department 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

A  simple  model  of  a  magnetostrictor  as  an  electro¬ 
mechanical  energy  conversion  device  is  presented  and  examined 
to  delineate  those  materials  variables  of  prime  importance 
in  magnetostrictive  transducers.  The  efficiency  of  the  model 
presented  is  examined  and  it  is  concluded  that  design  goals 
for  magnetostrictive  transducer  materials  should  be  to  obtain 
materials  lacking  porosity  and  internal  residual  stresses  and  ' 

having  a  large  magnetostrictive  strain  and  a  small  magnetic 
anisotropy . 

The  magnetostrictive  character  of  iron  rich  Fe-Co  and 
Fe-Al  alloys  is  presented  and  the  possibility  of  optimizing 
magnetic  anisotropy  and  magnetostrictive  strain  in  an  iron 
rich  Fe-Co-Al  alloy  is  discussed. 
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0.‘T  THE  FLOW  FIELH  AROUMb  A  MARI’TT?  nnCT.FO  PROPELLER 

by 

J.  A.  Fillo 
n.  F.  Cromack 

Mechanical  and  Aerospace  Engineering 
University  of  Massachusetts 
Amehurst,  Massachusetts 

This  report  describes  the  research  currently  being  con¬ 
ducted  on  marine  ducted  propellers  at  the  University  of  Massa¬ 
chusetts  under  Project  THEMIS.  A  brief  summary  of  the  progress 
is  included. 

The  experimental  test  facility  that  has  been  designed  and 
is  being  constructed,  along  with  some  of  the  instrumentation 
and  tests  to  be  conducted  are  described.  A  summary  of  experi¬ 
mental  and  analytical  methods  used  by  other  investigators  for 
ducted  propellers  is  presented.  Methods  of  improving  and  ex¬ 
tending  these  analyses  are  also  discussed.  A  theoretical 
approach  to  the  analysis  of  ducted  propellers  appears  to  be 
more  beneficial  than  an  experimental  parametric  study  because 
of  the  enormity  of  doing  the  latter.  Thus  the  theoretical 
results  are  eventually  to  be  compared  to  selected  experimental 
conditions . 
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PHYSICAL  AND  NUMERICAL  F.XPERTMF.NTC 
IN  HYDRODYNAMIC  WAKE  FLOWS 

by 

Thomas  J.  Mueller 
Aerospace  &  Mechanical 
Engineering  Department 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

Experimental  and  numerical  studies  were  performed  to 
determine  the  structure  of  the  viscous  wake  flow  behind  a 
backward  facing  step.  The  experimental  program  was  performed 
in  a  hydrodynamic  towing  facility  using  two  backstep  models. 
The  first  model  had  a  distance  from  nose  to  backstep  of  10 
inches,  a  step  height  of  0.50  inches  and  a  step  aspect  ratio 
of  6.  The  second  model  was  exactly  the  same  as  the  first 
model  except  that  it  had  a  step  aspect  ratio  of  9.  Both 
models  measured  9  inches  from  the  bottom  of  the  step  to  the 
end  of  the  splitter  plate. 

The  backstep  model  was  submerged  in  the  hydrodynamic 
twoing  facility,  and  rigidly  attached  to  a  motor  driven  camera 
cart  above  the  tank.  Aluminum  powder  was  suspended  in  the 
water  so  that  streak-time  photographs  could  be  taken.  From 
a  large  number  of  streak-time  photographs  and  using  an  optical 
comparator,  it  was  possible  to  determine  the  velocity  distri¬ 
bution  in  the  boundary  layer  upstream  of  the  step  as  a 
function  of  Reynolds  number.  These  velocity  distributions 
formed  part  of  the  input  for  the  numerical  studies.  The 
streak- time  photographs  also  produced  data  for  the  length  of 
the  separated  flow  region  to  reattachment  as  a  function  of 
step  Reynolds  number. 
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The  variation  of  measured  boundary  layer  thickness  with 
Reynolds  number  based  on  length  from  the  leading  edge  of  the 
model  wa s  compared  with  the  values  obtained  by  Blasius,  i.e., 
5/x^R  As  R  decreases,  the  Blasius  result 

SX  G  X 

that  6  increases  for  a  fixed  X.  The  measured  <5  approximately 
agreed  with  the  Blasius  result  from  R  of  7000  to  about  2000. 

G  X 

As  R  decreased  below  2000  to  the  lowest  value  obtained  of 

ex 

400,  however,  the  measured,  value  of  6  decreased  while  the 
Blasius  6  increased  sharply. 

For  each  Reynolds  number,  the  measured  velocity  distri¬ 
bution  was  inserted  as  part  of  the  inflow  boundary  condition 
for  the  numerical  study  of  planar  incompressible  laminar 
flow  over  a  backstep.  This  numerical  method  used  the  ap¬ 
propriate  time-dependent  finite  difference  equations  with 
physically  meaningful  boundary  conditions.  A  direct  comparison 
of  the  numerical  results  with  the  corresponding  experiments 
was  made  for  the  variation  of  length  of  the  separated  flow 
with  step  Reynolds  number.  As  the  step  Reynolds  number  in¬ 
creased  from  27  to  about  200,  the  dimensionless  length  cf  the 
separated  region  increased  for  both  the  numerical  and  experi¬ 
mental  results.  Furthermore,  this  comparison  indicated  a  very 
good  quantitative  agreement. 
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TRANSMISSION  OF  SOUND  THROUGH  WIRE  SCREENS  IN  STEADY  FT.OW 

by 

G.  A.  Russell 
S.  N.  B.  Murthy 

Mechanical  and  Aerospace  Engineering 
University  of  Massachusetts 
Amherst,  Massachusetts 

The  transmission  of  sound  across  a  wire  grid  located  cross- 
sectionally  in  a  small  tube  carrying  a  low  velocity  flow  is 
studied.  A  brief  history  of  this  and  related  problems  is  pre¬ 
sented  together  with  a  technical  justification  for  the  oresent 
research. 

The  physical  arrangement  of  the  wire  grid  and  tube  is  pre¬ 
sented  and  the  assumptions  necessary  for  writing  linearized 
perturbation  equations  are  discussed.  The  equations  are  given 
and  a  derivation  of  a  subsequent  equation  involving  the  steady 
flow  velocity,  screen  constant,  acoustic  frequency,  and  sound 
transmission  coefficient  is  outlined.  The  interpretation  of 
this  equation  in  terms  of  the  limiting  case  of  no  flow  is  dis¬ 
cussed.  Numerical  solutions  of  the  equation  are  presented 
and  briefly  interpreted  in  terms  of  air  flow.  An  extension 
to  the  case  of  other  flowing  fluids  is  also  outlined. 

Experimental  results  obtained  to  date  are  presented  and 
some  of  the  measurement  problems  are  discussed.  Plans  for 
future  analytical  and  experimental  studies  and  possible  appli¬ 
cations  for  deep  submersible  vehicles  are  given. 
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LAMINAR  FREE  CONVECTION 
FROM  A  NON- ISOTHERMAL  VERTICAL  PLATE 
IN  A  TET1PERATURE-STRATIFIEL)  ENVIRONMENT 

by 

K.  T.  Yang 

Aerospace  &  Mechanical 
Engineering  Department 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

The  temperature-stratified  ocean  environment  may  signi¬ 
ficantly  influence  the  rate  of  heat  transfer  from  a  heated 
surface  and  thus  the  performance  of  a  practical  heat-rejection 
system.  Such  temperature  stratification  also  exists  in  under¬ 
water  habitats.  In  this  study  the  phenomenon  of  laminar 
free  convection  along  a  non-isothermal  vertical  plate  in  a 
temperature  or  density-stratified  environment  is  investigated 
both  analytically  and  experimentally. 

Analytically ,  a  similarity  solution  in  terms  of  three 
physical  parameters  has  been  obtained.  Two  of  these  parameters 
refer  to  the  wall  and  ambient  temperature  variations,  while 
the  third  is  the  Prandtl  number  of  the  fluid  in  the  environ¬ 
ment.  Extensive  numerical  solutions  to  the  resulting  ordinary 
differential  equations  have  been  carried  out  for  wide  ranges 
of  these  parameters  which  are  of  physical  importance.  To 
help  determine  the  complete  effect  of  Prandtl  number,  efforts 
have  been  directed  toward  obtaining,  by  asymptotic  analyses, 
limiting  behaviors  for  Prandtl  numbers  approaching  zero  or 
infinity.  All  numerical  results  are  now  complete,  and  have 
provided  the  basis  for  discussing  the  physical  effects  of  the 
parameters  or  the  heat-transfer  characteristics  at  the  plate, 
as  well  as  for  demonstrating  the  fact  that  approximate 
solutions  based  on  quasi-uniform  environmental  temperatures 
are  not  adequate  to  predict  the  effects  of  ambient  temperature 
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variations  and  the  analysis  must  include  temperature  stra¬ 
tification  in  the  environment.  Furthermore,  several  int^re»L 
ing,  but  unexpected  behaviors  dealing  with  the  interaction  of 
parameters  have  been  noted,  and  have  been  explained  in  terms 
of  the  physical  processes. 

In  the  experimental  study,  an  environmental  tank,  60"  in 
height  and  30"  in  diameter,  has  been  designed,  fabricated,  and 
will  become  operational  shortly.  An  6"  Mach-Zehnder  inter¬ 
ferometer  purchased  by  matching  University  and  National  Science 
Foundation  funds,  has  been  ordered  and  this  interferometer 
will  be  used  to  map  detailed  temperature  or  density  field 
surrounding  the  heated  test  vertical  plate  mounted  in  the 
environmental  tank  through  8"  optical  flat  windows.  The  design 
of  the  tank  stand,  test  vertical  plate  and  instrumentation  is 
still  progressing. 
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THE  APPLICATION  OF  DYNAMIC  RELAXATION 
TO  THE 

FINITE  ELEMENT  METHOD 
OF 

STRUCTURAL  ANALYSIS 
by 

S.  Kelsey 

Civil  Engineering  Department 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

Dynamic  Relaxation,  as  developed  by  Day,  Otter,  Cassell 
et  al ,  is  a  matrix  iterative  method  for  the  solution  of 
simultaneous  linear  equations.  The  application  of  the  tech¬ 
nique  to  problems  of  structural  analysis  under  static  stress 
conditions,  as  used  by  these  investigators,  has  been  almost 
exclusively  associated  with  the  finite  difference  formulation 
in  space  of  both  the  equations  of  motion  and  the  constitutive 
relationships . 

In  this  report  the  authors  discuss  the  constraints  which 
the  earlier  investigators  have  placed  on  the  method  of  Dynamic 
Relaxation  by  restricting  themselves  to  the  finite  difference 
approach,  particularly  in  relation  to  the  static  stress 
analysis  of  plates  and  shells  having  arbitrary  geometric 
discontinuities  and  boundaries.  To  overcome  this  problem, 
an  alternative  approach  using  finite  elements  in  space,  in 
lieu  of  finite  differences,  is  proposed  and  developed. 

The  mathematical  basis  for  the  application  of  Dynamic 
Relaxation  to  the  solution  of  the  equations  of  motion,  and 
the  constitutive  equations  is  given  for  the  finite  element 
approach  in  space.  The  convergence  of  the  method  is  analyzed 
by  transforming  the  process  into  a  standard  eigenvalue  problem 
for  the  error  vector,  and  its  dependence  on  the  conditioning 
of  the  equations  is  demonstrated.  Values  of  the  parameters 
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used  in  the  iteration  to  optimize  convergence  are  derived.  The 
theory  for  the  optimum  asymptotic  convergence  of  Dynamic 
Relaxation  is  shown  to  be  applicable  to  a  much  wider  range 
of  problems  than  is  the  method  of  Successive  Overrelaxation. 

For  the  special  case  of  a  tridiagonal  form  of  the  coefficient 
matrix,  a  comparison  of  the  optimum  convergence  of  Dynamic 
Relaxation  is  made  with  that  of  the  other  basic  iterative 
methods,  namely  Point  Jacobi,  Gauss-Seidel,  and  Successive 
Overrelaxation.  In  this  particular  instance,  the  convergence 
of  Successive  Overrelaxation  is  shown  to  be  better  than  that 
of  Dynamic  Relaxation. 

A  strategy  for  approaching  the  optimum  convergence  of 
Dynamic  Relaxation  is  developed  and  demonstrated  by  the  appli¬ 
cation  to  plane  stress  problems  which  include  the  static  stress 
analysis  of  flat  plates  with  circular,  elliptical  and  filleted 
square  holes. 
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THE  NONLINEAR  ANALYSIS  OF  PLATES  &  SHELLS 

by 

Dr,  M.  M.  Miller 
Dr.  A.  Chajes 

Department  of  Civil  Engineering 
University  cf  Massachusetts 
Amherst/  Massachusetts 

Using  the  finite  difference  technique,  a  method  for  non¬ 
linear  shell  analysis  is  being  perfected.  Newton's  method  is 
used  to  linearize  the  nonlinear  terms  in  the  differential 
equation.  A  load-deflection  curve  is  obtained  by  increasing 
the  load  in  small  increments.  For  each  increment  the  variable 
coefficients  that  arise  from  Newton's  expansion  are  obtained 
by  an  iteration  procedure.  Consequently,  a  large  set  of 
simultaneous  linear  algebraic  equations  is  solved  repeatedly. 
If  the  behavior  of  the  shell  is  described  by  second  order 
differential  equations  one  obtains  a  tri-diagonal  block  system 
after  substitution  of  the  difference  ratios.  The  solution  is 
obtained  by  Gaussian  elimination  where  the  individual  elements 
are  blocks  each  consisting  of  several  unknowns. 

Using  a  shell  with  double  curvature  the  above  algorithm 
has  been  tested  and  found  to  be  satisfactory. 

The  variables  in  the  second  order  differential  equations 
are  inplane  and  transverse  displacements.  This  formulation 
is  useful  only  if  the  boundary  conditions  are  given  in  terms 
of  displacements.  If  the  boundary  conditions  are  given  in 
terms  of  forces,  a  formulation  consisting  of  two  second  order 
equations  becomes  more  convenient.  Using  a  plate  loaded  both 
axially  and  transversely,  various  boundary  conditions  were 
investigated.  Both  the  second  order  equations  and  the  fourth 
order  equations  gave  satisfactory  results  when  used  in  cases 
with  appropriate  boundary  conditions. 


For  a  shell  subject  to  instability  more  than  one  equi¬ 
librium  configuration  may  exist  in  the  large  deflection  domain. 
It  is  possible  that  a  shell  may  "snap  through"  from  one  of 
these  to  another.  This  phenomenon  exists  also  in  plates  and 
has  been  studied  extensively  for  that  case. 

The  methods  described  above  are  now  being  applied  to  a 
cylindrical  shell  under  external  pressure  and  axial  compression. 
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INELASTIC  AXISYMMETRIC  BUCKLING 
OF  RING  STIFFENED  CYLINDRICAL  SHELLS 

by 

L.  H.  N.  Lee 

Department  of  Engineering  Science 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

In  this  paper,  the  effects  of  prebuckling  deformation 
caused  by  the  ring  constraints  on  the  inelastic  axisymmetric 
buckling  of  ring  stiffened  cylindrical  shells  under  external 
pressure  are  studied.  By  Drucker's  postulate  of  positive 
work  in  plastic  deformation,  a  modified  incremental  theory  of 
plasticity  on  the  stress-strain  relationships  expressed  in 
terms  of  Kirchhoff  stress  and  Green  strain  rates  is  developed. 
Furthermore,  a  variational  principle  in  Lagrangian  descriptions 
and  for  quasi-static  problems  of  finite  plasticity  is  developed 
and  employed  to  solve  the  problem. 

It  is  shown  that  the  modified  constitutive  relation¬ 
ships  satisfy  the  coordinate  and  spatial  invariance  and  energy 
requirements.  Using  the  constitutive  relationships,  the 
existence  of  an  extremum  principle  for  a  material  having  a 
sufficiently  great  rate  hardening  is  shown  and  a  criterion 
for  »he  stability  of  a  body  under  dead  load  is  established. 

The  variational  principle  in  conjunction  with  the  incre¬ 
mental  Rayleigh-Ritz  technique  is  used  to  determine  the  de¬ 
formation  process  of  an  idealized  cylindrical  shell  composed 
of  four  thin  load-carrying  sheets  made  of  a  general  work 
hardening  material.  A  comparison  of  the  theoretical  pre¬ 
dictions  obtained  numerically  by  an  Univac  1107  digital 
computer  with  available  experimental  results  is  made.  The 
comparison  is  favorable. 


LIGHT  SCATTERING  AS  A  METHOD  OF  DETERMINING  THE 
ACOUSTICAL  AND  TRANSPORT  PROPERTIES  OF  FLUIDS 

by 

A.  A.  Ilonkewicz 

Department  of  Engineering  Science 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

A  high  sensitivity,  high  resolution  system  for  the  study 
of  light  scattered  from  thermal  fluctuations  in  gases  has  been 
constructed.  A  Bendix  Channeltron  photomultiplier  tube  with 
a  dark  noise  of  one-tenth  of  a  photoelectron/sec  is  the  basis 
of  an  experimental  setup  which  includes  a  Fabry-Perot  inter¬ 
ferometer  having  a  band-width  of  10  MHz  with  a  piezo-electric 
scan  and  single  mode  He-Ne  laser.  Resolution  of  10  MHz  is, 
therefore,  attainable  and  signal  levels  as  low  as  four-tenths 
of  a  photoelectron/sec  yield  reliable  data. 

With  this  apparatus  the  density-density  correlation 
function,  G(r,t) ,  of  several  gases  has  been  investigated  for 
a  range  of  pressures.  Whether  the  form  of  G(r,t)  is  derived 
from  a  kinetic  or  hydrodynamic  view  depends  on  the  value  of 
the  Knudsen  number,  i.e.  the  ratio  of  the  acoustic  wavelength 
to  the  mean  free  path.  Previous  studies  have  concentrated 
on  varying  the  angle  of  scattering  (Changing  the  acoustic 
wavelength)  while  this  study  varies  the  gas  pressure  itself 
(changing  the  mean  free  path)  allowing  a  more  detailed  in¬ 
vestigation  of  the  transition  from  the  kinetic  to  the  hydro- 
dynamic  view. 
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DYNAMIC  PROGRAMMING 
AND  THE  DISCRETE  MAXIMUM  PRINCIPLE 

by 

R.  J.  Leake 

Department  of  Electrical  Engineering 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

As  with  continuous  dynamical  system  optimization,  the 
basic  approaches  to  optimal  control  of  discrete  time  systems 
have  their  origin  in  the  variational  theory  of  mechanics; 
Dynamic  Programming  stemming  from  the  Hamilton-Jacobi  Theory 
and  the  Discrete  Minimum  Principle  based  on  the  Hamilton 
Canonical  Equations. 

Necessary  and  sufficient  conditions  for  optimality  of 
discrete  time  systems  over  arbitrary  sets  are  given,  followed 
by  some  interesting  observations  concerning  global  successive 
approximations.  Examples  are  given,  and  connections  between 
the  minimum  principle  and  dynamic  programming  are  established. 
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MICROSTRUCTURE  EFFECTS 
IN  THE  ANALYSIS  OF  COMPOSITES 

by 

T.  Ariman 

Department  of  Engineering  Science 
University  of  Notre  Dame 
Notre  Dame,  Indiana 

In  this  report  the  theory  of  micropolar  continua  has  been 
used  to  analyze  the  stresses  and  displacements  within  the 
separate  components  of  a  fibrous  composite  material.  The 
separate  constituent  materials  are  assumed  to  be  linearly 
micropolar  elastic,  isotropic,  and  homogeneous.  A  model 
chosen  consists  of  parallel  fine  fibers,  embedded  in  a  matrix 
material.  The  problem  is  formulated  as  one  of  plane  micro- 
polar  elasticity  and  solved  for  some  special  cases.  Various 
numerical  results  are  presented  and  compared  to  those  of  the 
classical  and  indeterministic  couple  stress  theories. 
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